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^ Abstract 

A Quantum Field Theory formulation of Bose- Einstein Correlations is given. 
^ It contains as a special case the classical current approach. It is shown that 

^ the particle-antiparticle correlations are a general feature of Bose-Einstein 

^~~^ Correlations and not an artifact of a certain approximation. Relations are 

O^ provided through which the quantum corrections to the classical current 
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approach can be determined experimentally. 
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Bose- Einstein correlations (BEC) are an important theoretical and experimental tool 
in particle and nuclear physics in a wide range of energies because one can obtain from 
BEC information about the space-time structure of the source as well as about the dy- 
namics of the reaction mechanism. 
There have been two main theoretical approaches to BEC: 

1. The wave function approach [1]. This approach being based on first quantisation 
has serious drawbacks some of which are listed below: 

(a) The treatment of inclusive reactions which constitute the bulk of experimen- 
tally accessible reactions is quite cumbersome, if not practically impossible for 
high multiplicities. 

(b) Higher order correlations can hardly be calculated with this method. 

(c) A correct treatment of quantum statistical coherence is not possible. 

(d) Quantum statistical particle- antiparticle correlations are absent. 

2. The classical current formalism [1, 2], which has none of the drawbacks mentioned 
above but still uses the approximation that the currents (but not the fields) are 
classical. The Wigner function approach [3] and the method of local equilibrium 
statistical operator [4] are particular cases of this current ensemble approach. 

Recently using a particular model [5] for heavy-ion reactions it has been conjectured that 
measurements of a classical radius through BEC are strongly affected by non-classical 
off- shell effects. 

In this paper we develop a quantum-field approach to BEC in the real time formula- 
tion (Keldysh technique) [6]. Among the results obtained here we mention an estimate 
accessible to experimental measurements of the quantum corrections to the classical cur- 
rent formalism and the proof that quantum-statistical particle- antiparticle correlations are 
not an artifact of the classical current formalism but a quite general property of particle 
physics. 

We start with the inclusive spectra which are usually defined in terms of the density 
matrix pj of the final state 

J2^m J. J. 

P2(fci, ^2) = kXirr-irr- = Tr{pfJik2)Jik,)aik,)aik2)} = (2) 



(PkicPk 
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[ki^k2) ; ^jWfnf{ki)nf{k2) =<n{ki)n{k2) > 

[ki = k2 = k) ; E/ WfUfik) ■ {nf{k) - 1) =< n{k) ■ {n{k) - 1) > 



(where wj denotes the probabihty to obtain the final state |/> from the initial state and 
n is the multiplicity) and express them in terms of the initial state density matrix pi and 
the S-matrix using the identity pj = S piS"^ 

Pi(fc) = Trp,sK^k)a{k)S (3) 

P2(fci, fc2)=fci°fc° 3. „. = Tr{p,sKHk2)Jik,)aik,)aik2)S} (4) 

To proceed further we assume [7] that in our interferometry experiment we registrate 
particles with quantum numbers not contained in the initial state 

ap, = p,a'^ = (5) 

This condition is in general satisfied if one omits in the final state particles with zero 
transverse momentum. We can transpose the annihilation operators a on the r.h.s. and 
the creation operators aT on the l.h.s. side of Tr expression commuting repeatedly a with 
S, aT with S'T and using (5). This results in the following expressions for the single 

P^ik) = TrL.^^-^] (6) 

^ ' Y Sa{k)sJ{k)j ^ ' 

and the double inclusive cross-section 

''6a{k,)6a{k,)-sJik,)SJ{k,)] ^^^ 

In (6) and (7) we used the property [a{k), S] = SS/Sa'^{k), where the variational derivative 
is taken with respect to the corresponding symbol of the operator aT (cf. [8]). 

To be specific we consider in the following the direct production of pions Tr^"'"^ Tr^"^ 7r° 
(the method can be easily generalized for pions coming from resonances). 
We assume that the particle production is governed by the interaction lagrangian: 

iv™*(x) = j(+)(x)7r(-)(x) + j(-)(x)7r(+)(x) + J°(x)7r°(x) (8) 

which is the simplest ansatz for an isospin triplet. The currents J^~^\ J ■, J^ are 
operators not depending on the tt^^' and 7r° fields ([J, tt] = 0). For example J can be due 
to quark-fields {J^'{x) = u^sd , J° = (u^su — d^^d)/ y2) or any other fields different 
from the pion fields. In this way we generalise the classical current formalism. 
Let a(aT) , h{lr) , c(cT) stand for annihilation (creation) operators of tt^' , tt^~^' , tt^^' 
respectively. The same procedure that was used to derive (6), (7) leads now to following 
expressions for inclusive spectra: 



p(o.W.«.°^ = r.L.i^^| (11, 

' ^"''"'^-'''^d-k.d-k, - ''Y'6a{k,)6a{k,) 6a^{k,)6a^{k,)j ^ ^ 
' ^ ^' ^^- ^ ^^^M^fc2 ~ YSc{k,)6c{k,)- ScHk,)6cHk,)j ^ ^ 

' ^"''"'^-''''d^k.d^k, - ''Y6a{k,)6b{k,) sJ{k,)6bHk,)j ^ ^ 

Formulae (9)-(ll) represent single inclusive cross-sections of tt^") , tt^"'") and 7r°. Double 
inclusive spectra for 7r~7r~ = tt^"'") tt^"'") and ir^^' ir^^' pairs are described by (12) and (13) 
respectively. Formula (14) stands for the tt^") tt^"'") double inclusive cross-section, which 
reflects Bose-Finstein Correlations between a particle and its antiparticle [2]. 
Next we eliminate the field operators in favour of the current operators J by means of 
the S'-matrix represented through the interaction lagrangian 

5" = Texpjz /j^xiv™*(x)| (15) 

and calculate the derivatives: 

6S T{j(-\k)S} 6S^ f{sU^+\-k)} 



sJ{k) 


V2Py 


ss 


T{ji+Hk)S} 


Sb^k) 


V2Py 


ss 


T{j(^\k)S} 



Sa{k) ^Fy 

,55't _ f{5'tj(-)(-fc)} 

Sb{k) ~ V2W 

SS"^ t{5'tj(o)(-A;)} 



(16) 



6c^{k) V2m/ ' Sc{k) y/Wv 

Here T and T are the chronological and antichronological time ordering operators re- 
spectively. No difficulties arise for the second and higher order derivatives because the 
currents J do not depend on a, aT,6, fcT^ c, cT. 
The single and double inclusive spectra read now: 

p[-\k) = Tr L,f[sU^+\-k)]T[j^-\k)S] I (17) 



Pl+\k) = Tr L,f[sU^-\-k)]T[j^+\k)S] I (18) 

pf)(fc) = TrL.f [5'tj(°)(-A;)]T[j(°)(A;)5']| (19) 

Pi--)(fci,fc2) = rr{/>,f[^tj(+)(_^^)j(+)(_^2)]T[j(-)(^i)'/^"H^2)^]} (20) 
Pr^(fci,fc2) = rr{^.f[^tj(o)(_^^)j(o)(_^^)]^[j(o)(^^y(o)(^^^^]| (21) 

Pi-+)(fci,fc2) = rr{/>,f[^tj(+)(_^^)j(-)(_^2)]T[j(-)(^i)'/^+H^2)^]} (22) 

We have thus obtained the expressions for the physical quantities only in terms of the 
currents and S'-matrix. To eliminate the S'-matrix we rewrite (17) - (22) in terms of 
Heisenberg operators. To do so we insert the identity 1 = SS^ between the antichrono- 
logical T and chronological T parts of expressions (17) - (22) and use the well known 
relation between Heisenberg and Interaction representation operators: 

T [AhBhCh ■ ■ ■] = S^T [A,ntBintC\nt ■ ■ ■ S] 

and (23) 



tJ A-in,t-Din,tLjin.t ■ ■ ■ 



s 



T[AhBhCh...] = T 
This procedure results in the following formulae: 

Pi-\k) = Tr {pA^\-k)jii\k)} (24) 

pW(fc) = Tr {paii\-k)ji;\k)} (25) 

Pi'\k) = Tr\paS\-k)jS\k)} (26) 



Pt-\k, , k,) = Tr{p,f [j(+)(-^,)j(+)(-^2)] T [Ji,-\k,)ji,-\k,)]} (27) 

Pi''\k, , k,) = Tr {p,f [./^\-k,)./^\-k,)] T [jS\h)jS\k2)] } (28) 

Pt^\k, , k,) = Tr{p,f [j(+)(-^,)jH(-^2)] T [Jij-\k,)ji,^\k,)]} (29) 

Similar expressions can be derived of course for higher order inclusive cross-sections. 
From now on we shall omit the index H assuming that all operators are in the Heisenberg 
representation. Eqs. (24)-(29) constitute the Quantum Field theoretical formulation of 
the inclusive spectra in terms of quantum currents and the density matrix of the initial 
state. 



For the particular case of recoilless emission the quantum currents J can be replaced 
by classical currents [9] which simplifies to a great extent the expressions above. Thus 
e.g. for TT^~^' TT^' spectra one obtains as in [2] 

Pi;+)(fci, ^2) =[dJ^-^DJ^+^P{J^-\J^+^) X (30) 

where the trace over the density matrix pi has been replaced by an integral over the 
classical current configurations. 

In Quantum Field Theory there is a very elegant formalism to deal with antichronolog- 
ically ordered operator products on the same footing as with chronological ones, namely 
the Keldysh technique [6] well known also under the name Real Time Formalism (RTF) 
The essence of this formalism consists in the following: 

Fvery dynamical operator A is mapped formally in two operators denoted here as A(^ , Aq. 
These additional degrees of freedom refiect what kind of ordering is used (Fx.: Tc[AqBq] = 
T[AB] ; Z[A^B^] = f[AB] ; Z[A^Be] = AB ; Z[AeB^] = BA). In RTF the double- 
inclusive spectra (27)-(29) are just ^^ contour-time^^ (Tc) ordered Green-function of the 
corresponding operators: 

-k,)jl^\-k,)jl^\k,)jl^\k,)]> (31) 

-k,)J§\-k,)jl^\k,)jl^\k,)}> (32) 

-k,)jl^\-k,)jl^\k,)jl^\k,)]> (33) 

Above we have introduced the notation <...>= Tr{pi . . .}. Formally these equations 
are quite similar to those obtained in the classical current formalism (compare e.g. Fq. 
(33) with Fq. (30)); the difference resides, of course, in the fact that in the quantum-field 
formulation the currents are operators and they appear in a particular order. 

It is useful to rewrite the full Green-functions (denoted by <...>) in terms of their 
connected parts (denoted here as ^ . . . ^). For simplicity we assume from now on that 
the system is completely chaotic (< J >= 0). Then, the double-inclusive spectra (31)-(33) 
written in terms of connected Green-functions take the following form (for their graphic 
representation (cf. Fig.l): 

Pt~\k^,k,) = <T,{4+)(-^i)4^V^2)4"H^i)4"^(^2)}> + 

+ <T, {4+V^i)4"^(^i)}><'^c {4+V^2)4"^(A:2)}> 

+ <T, {4+V^i)4"^(^2)}><T, [jl^+\-k,)j!-\k^)}:> (34) 



i--\k^,k2)-- 


-<z{f^' 


pr\k„k,)- 


-<z{jS\ 


i-+\k,,k2)-- 


-<Tjjt^ 



p^''\krM) = <T,{j§\-h)j§\-k,)ji;\h)ji;\k,)}:> 

+ <T, {4°)(-^i)4°^(^i)}><'^c { jr(-^2)4°\^2)}> 
+ <T, {4°)(-^i)4°)(^2)}><T, { j(°)(-^2)4°^(^i)}> 

+ <T, { j(°)(-^i) j(°)(-^2)}><T, {4°)(^2)4°\^i)}> (35) 

Pt^\k^,k2) = <:z{ji^\-h)jl-\-h)4-\h)4+\h)}:> 

+ <T, {4+V^i)4"^(^i)}><'^c {4" V^2)4^^(A:2)}> 

+ <T, [,4+\-k,)jl-\-k2)}:><:Z {4~\k^)4+\k2)}:> (36) 

The expressions (34)- (36) are exact but quite complicated. 

An important simplification can be achieved by assuming a Gaussian form for the 
density matrix [10]. This form follows from the Central Limit Theorem which is appli- 
cable when the number of independent sources becomes large. The assumption that the 
number of sources is large is quite natural for high energy heavy-ion reactions and it is 
an accepted ingredient in many phenomenological approaches like e.g. the string model 
[11, 12], hydrodynamical [13] and thermodynamical models, etc. For a gaussian density 
matrix the 4-point connected correlator vanishes. Then eqs. (24) - (26), (34)-(36) become 

P[^\k) = F^{k,h) (37) 

p[-\k) = P[-^\k) = F'''{k,k) (38) 

Pt~\k^,k2) = Pt\k^)Pt\k2) + \F^''{k^,k2)\' (39) 

P^''\k,,k2) = Pi'\k,)Pi'\k2) + \F-^ik,,k2)\' + \F-^ik,,k2)\' (40) 

Pt^\k,,k2) = Pi-\k,)Pi+\k2) + \F^^i-k,,k2)\' (41) 

where the functions F are defined for charged particles (upper index ch) and for neutral 
ones (upper index n) as follows: 

F^'^ik^^k^) = <Tc {4^H-^i)4"H^2)}> = <.J^^\-k^)J^-\k2)'> (42) 

F'\k^,k,) = <Tc{jl^\-k^)4~\k2)]'> = <t{j(+)(-A;i)j(-)(A;2)}> (43) 

F-ik,,k2) = <Z { Jr(-^i)4°^(^2)}> = <j(°)(-^i) J(°)(^2)> (44) 



F-{h,h) = <T,{j(°)(-^i)4°^(^2)}> = <:f {j^'H-ki)J^'Hk2)}:> (45) 

In contrast to the classical current approach [2] which deals with only one type of two- 
current correlator we have here two different kinds of two-current correlators depending 
on their ordering prescriptions. The two-particle correlations are given by 

Ct^\k,,k,)^^S^/^^%^ = 1 + \d'\h.h)\^ (46) 



Ct~\k,,k,)= :_l '-y::'' ^ = 1 + \d^'{h,h)\' (47) 



Pt\k,)pi^\k,) 


Pt~\k,,k,) 


~ Pt\k,)Pt\k,) 


-f2 I'^U '^■2) 



^^°°^^^^'^^^^i#0#(M ^ ^ + \d''iki.k,)\' + \d-{h,h)\' (48) 

where the functions d are defined as follows: 

^"iM,k,). , ^"'^-^'' ; ^-iM,k,i. , ^"'^■■^-' ,49) 

^F"'{k„h)F"'{k„k2) ^F*(ki,k,)F*(k2,k2) 

rik,M^ , ^"'^■'^'' ; rjk^M^ , ^"'^■'^'' (50) 

^F'{k„ki)F'{k„k2) ^F'{k„ki)F'{k„k2) 

The "surprising" effects found in [2] i.e. the 

• presence of particle-antiparticle Bose- Einstein type correlations and 

• a new term in Bose- Einstein correlation function for neutral particles 

are reobtained, but under a more general form which contains also the quantum correc- 
tions. These equations also prove that the effects derived in [2] are not an artifact of the 
classical current formalism but have general validity. 

To study the deviations of the above equations from the classical current approach we 
compare the two types of two-current correlators (42) and (43). This can be done more 
easily in coordinate representation {F'^^{xi,X2) = ^J® {xi)Jq {x2)^ ', F'^^{xi,X2) = 
«t{jW(x04-^(x2)}») : 

R'\x^,X2) = F'\x^,X2)-F'\x^,X2) = Q{x^,-xl)<:[J^+\x^)J^-\x2)]:> (51) 

The i?^'' -function is just the retarded current commutator and can in principle be cal- 
culated once the currents and fields entering (8) are given (cf. the end of the paper). 
In principle, one can determine R'^^ experimentally by comparing 7r(~V(~) and Tr^"^^"'") 
data and using Eqs. (38), (46), (47), (49), (51). In this way one can test to what extent 
the system can be treated quasi-classically. On the other hand the i?-function can be 



calculated for a particular model. For any classical current theory i? = and therefore 
d{ki,k2) = d{h,-k2) [14]. 

Similarly for the neutral particles one has 

i?'^(xi,X2) = F'^(xi,X2)-F'^(xi,X2) = 0(x'|'-x°)<[j(°)(xi),j(°)(x2)]> (52) 

which is the retarded commutator of the current J^^>. 

The structure of formula (48) is typical for the correlation function of any neutral boson. 
For instance, if one measures photon-pairs with a given polarisation, formula (48) is appli- 
cable if one defines J^^'{k) = e^{k)J^{k). In general polarisation effects for particles with 
non-zero integer spin can play an important role in (40) [15]. But apart of polarisation 
effects the kinematical and dynamical information about the collision is contained in the 
two functions F"^ and F"^ . 

Since \d'^{k^k)\ = 1 and \d"'{ki,k2)\ > the intercept of the correlation function C2 {k^k) > 
2 (we refer here to a completely chaotic source, without coherence). The exact value of 
the intercept is model dependent. As mentioned above, in the classical current formalism 
R = 0, d'^[ki^k2) = d'^[ — ki^k2) and the correlation function simplifies to 

C£°)(fci,fc2) = l + \d-{k^,k2)\' + \d-{-k^,k2)\' (53) 

which depends only on the function J". For this case the bounds for the correlation 
function have been established [16]^. 

To exemplify the above made considerations we present here a simple model where 
massless pions are produced from a static source consisting of other scalar particles in 
thermal equilibrium (p = e~^'^ where H is the Hamiltonian of the system and T is the 
temperature, cf. e.g. [18]). The interaction is given by the Lagrangian 

L^nt = gj{x)Tr{x) , with j{x) = ip{x)ip{x) (54) 

The calculations are performed in the one-loop approximation using the Keldysh technique 
[6]. 

The relevant current correlators can be expressed in a very simple form if one takes 



^In [17] the results of the classical current formalism were misquoted. Contrary to what was stated in 
[17] the intercept of the tt^tt^ correlation function in the classical current formalism strongly depends on 
both the mass m of the particle and the time-duration of freezeout r. The value 3 for the intercept can 
be obtained either in the limit of instantaneous freezeout or for massless particles at zero momentum. 
For the "static" source with exponential parametrisation one has C™^^ = 2 + exp(— rm) . 



into account the contour time ordering by placing (H — ) labels at the vertices: 



< jixi)jix2) > 



< j{x2)j{xi) > 



+ 



Xi 



Xi 



X2 



+ 



X2 



(< ^{XIMX2) >f 



v< ^(X2)'^(Xi) >] 



^^""''^ < ^ix,Mx2) >=\ I 7^e-^(— ^)27r,5(/)[0(po) + n{\p,\ 



(55) 



h J (27r) 

The i?-function which is the retarded current commutator (52) can be calculated even 
easier than the correlators 



[j{xi),]{x2)] = [<~p{xi)Lp{xi) , Lp{x2)<~p{x2)] 

= Cp{xi)[Lp{xi) , Lp{x2)\^{x2) + ^{x2)[Cp{xi) , Lp{x2)\^{xi) 

because the field commutator [',c(xi) , Cp{x2)\ does not depend on temperature: 

[^(^i), ^(^2)] = ^/^e-^(---)27rsgn(po)^(/) 
We omit the algebraic details and present just the results: 



R[xi^ X'l 



-\g'hTf{x,)Q{x',-x^)^-^f{x2) 
27r^ \ux\ 



1 1 r 

32^» - 1 ~ 27ra ^ 2. 



(56) 



(57) 



(58) 



< ]{xi)]{x2) >= f{xi)g^ 



h 



2^"^ x"^ — i€sgn[xo) 



+ 



(59) 



T ^r 

— / da sgn{a — {ux))S i{au'^ — X 
27r J ^ 



A* ^A'\2 



1 

+ 77 



^2naT/n _ I 27raT/n 2 



fix2 



where x = Xi — 2:2 and we define a = 2T\ux\/h using the 4- velocity u^ ( in rest frame 
M^ = (1, 0)). The function f{x) is introduced phenomenologically and plays the role of an 
effective vertex function containing the finite size effects of the source. It is the subject 
of Bose- Einstein correlations studies. 



It is easy to see from (58) and (59) that this model manifests quasi-classical behaviour 
(that is \R{xi,X2) l< j{xi)j{x2) >\ <^ 1) at least in two limiting cases: T — )• oo or /i — )• 
as it should be from general principles. 

Although this model can not be considered as realistic it illustrates how quantum effects 
do influence BEC. 

Finally one should stress a new feature of particle-antiparticle BEC which emerges 
from this paper. Besides the fact that like ordinary particle-particle BEC they can serve 
as tools in the determination of radii and lifetimes, and besides the fact that they reflect 
the existence of squeezed states, they play a particularly important role in the detection 
of quantum corrections to BEC. Their experimental measurement is a highly rewarding 
task. 

We would like to acknowledge fruitful comments by I.V.Andreev, C.Jeffries, H.Paul, 
M.Pliimer. This work was supported in part by the Gesellschaft fiir Schwerionenforschung, 
Darmstadt. 
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+ 

ki kt • 


!>• ki —ki 






h -h • 


D^ k2 -k2 



+ 




ki 



+ 



+ 



Fig. 1 Decomposition of the ir^ > ir^ > 4-point full Green-function. 



in terms of the connected parts. This diagram corresponds to P2 (^i? ^: 



(—), 



-h 



-h 



+ 




k\ —k\ 



ko —k- 



Q 



kx -kx 

+ 

k^ —ko 



+ 



k\ —k\ 

+ 

k^ —k^ 



+ 




kx -kx 

+ 

k^ —ko 



+ 



ki 



Fig. 2 Decomposition of the ir^^' ir^^' 4-point full Green-function 

in terms of the connected parts. This diagram corresponds to P2 (ki, k2). 



-kr 



+ 




ki_ -k-i 



ki -h 



+ 



+ 




+ 


+ - 






I 




— ^ k, -ko V i 



ki 



+ 



+ 



+ 



Fig. 3 Decomposition of the tt^ > tt^"'") 4-point full Green-function 



(-+), 



in terms of the connected parts. This diagram corresponds to P2 (^i? ^: 
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